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A grand canonical system of hard-core bosons in an optical lattice is consid- 
ered. The bosons can occupy randomly N equivalent states at each lattice 
site. The limit N ^ oo is solved exactly in terms of a saddle-point integra- 
tion, representing a weakly-interacting Bose gas. In the limit N ^ oo there 
is only a condensate if the fugacity of the Bose gas is larger than 1. Correc- 
tions in 1/N increase the total density of bosons but suppress the condensate. 
This indicates a depletion of the condensate due to increasing interaction at 
finite values of N. 

PACS numbers: 03.75. Fi, 05.30.Jp, 32.80.Pj, 67.40.-w 



1. INTRODUCTION 

Recently developed systems of trapped ultra-cold gases open a wide new 
field for the study of many-particle physics. This includes the investigation 
of Bose-Einstein condensation of interacting Bose particlescla, leading to 
macroscopic quantum states or "matter waves" . Among the most interestir 
experiments there is the trapping of bosons in so-called optical latticesQ'E 
that are created by intersecting Laser fields. Typical lattice constant of 
optical lattices are a; = 250...500nm. They exhibit unsual properties due 
to the interference of macroscopic quantum states. The interaction between 
the atoms is characterized by the ratio of the scattering length a and the 
typical particle distance as an^/^, where n is density of particles. In 

most of the experiments we have an^^^ < 10~^ such that these systems are 
diluted. Interaction plays only a weak role in these systems such that the 
Gross-Pitaevskii approach is sufficientfl'i. After the experimental discovery 
of Feshbach resonances in the trapped Bose gases it became possible to tune 
the inter-atomic interaction with an applied magnetic field over a wide range 
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of interactionsi'0'0. The scattering length for ^^Rb atoms is a ~ 5nm but 
near Feshbach resonances it is up to 5OOnm0. In this case one can study an 
optical lattice in which the lattice constant (determined by the wave length 
of the applied Laser field) is equal to the scattering length. Then at most 
one boson can be found in the minima of the optical lattice. 

From the theoretical point of view this field requires the study of models 
with stroriely interacting particles. The so-called Bose-Hubbard model has 
been usedoH to decribe an on-site interaction with coupling strength U . A 
mean-field calculation! has shown that the condensate vanishes in the strong- 
coupling limit U — > oo, indicating that a strong interaction among the bosons 
has a destructive effect on the condensate. On the other hand, from the 
experiments it is believed that the interaction can still be considered as 
hard-core, characterized only by the scattering length a but including many- 
particle collisions. 

In the following we will tune the interaction by going from the weakly 
interacting (dilute) to the strongly interacting hard-core Bose gas. This will 
be achieved by the assumption that each lattice site can accommodate N 
bosons with the same energyEl. Bosons can randomly change between these 
states and tunnel between neighboring lattice sites from/to any of these N 
bosonic states with the same rate. Interaction exists only between bosons 
at the same state. The regime ^ 1 describes a dilute Bose gas whereas 
= 1 corresponds to a hard-core Bose gas with only one state per lattice 
site. We will solve the limit N ^ oo within a saddle-point integration and 
apply a 1 /A^-expansion to control systematically the corrections to the dilute 
Bose gas towards a more strongly interacting bosons. 



2. BOSONS IN AN OPTICAL LATTICE 

In the first part of this section the Bose gas on a lattice will be discussed 
without interaction using a complex field cp. This field will be replaced in 
the second part by a hard-core field which has the same propagator as the 
field (j). Therefore, in contrast to the Bose-Hubbard model, the interaction 
is carried by the field and not directly by the Hamiltonian. A model of 
non-interacting bosons is considered on a d-dimensional hypercubic lattice 
with M sites. A boson can occupy statistically one of A^ degenerate states 
at each site, and tunneling between these N states at lattice site r to any 
of the states at site r' occurs with the rate J°'? . This is represented by the 
Hamiltonian 

1 ^ 

^ = -^ E T.J^y'^t'-'^r' (1) 

a,a'=l r,r' 
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with bosonic creation and annihilation operators The statistics 

of a system of bosons can be described by introducing a grand-canonical 
ensemble with fugacity (" at the inverse temperature (3. In the case of non- 
interacting bosons it is defined by the partition function!! 



(2) 



where ej are the eigenvalues of the Hamiltonian matrix 



TTOl,a 



_1 

2 r,r' 



Assuming that (3 has integer values, the factors in Eq. (^) can be expressed 
by the determinant of a /? x /3-matrix as 



where 



/O u 
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Wn 



det(l — Qwj), 
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\u 
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... oJ 



with u = e . 



(3) 



Then the partition function can also be written as the determinant of a 
/37V X /3A/'-matrix as 

^ " det(l - (w) ■ ^'^^ 

It) is obtained from expression (P) by using for u the TV x TV-matrix u = 
. It is convenient for the following to introduce space-time coordinates 
X = {t,r) with a "time" variable t = 1,2,...,/?. Moreover, we assume that 
the matrix H = — log u is chosen such that 



N 



(5) 



2.1. Random- Walk Expansion 

Using a complex field (j)x the inverse determinant in Eq. (Q) reads 

Z = / exp [ - ^ ^ ^^^^ + ^ E E ^x,-"^"*/-"'] n #x#"/vr. (6) 



a; «=1 



x,a;' a,o'=l 
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The field is subject to periodic boundary conditions in t 

rp+l,r = ^l,r- (7) 

The integrand of Z in the expression (^) can be formally expanded in powers 
of the matrix w as 

a,a'=l t=l r,r' a,a'=l l>0 

(8) 

The integration over the Bose field (px can be performed for each expansion 
term, leading to a random- walk expansion of Z with elements Cu^.r'-I^l^ 
From Z we can evaluate the density of bosons as 

In the case of the non-interacting Bose gas Eq. (^) we get immediately 

w'' contains I — 1 summations with respect to a = 1,...,N, contributing a 
factor A^'~^. Another factor A'^ comes from the trace. Altogether, this gives 
n oc A^~^, a consequence of the scaling in Eq. (||). It should be noticed that 
C is restricted here to < 1 if the eigenvalues w/N are < 1. This is a 
well-known artefact of non-interacting bosonsH. As we will see subsequently, 
> is not restricted in the interacting Bose gas and has a non-zero density 
for C > 1 even in the limit N ^ oo. 



2.2. The Hard-Core Bose Gas 

A hard-core condition can now be implemented in this model by assum- 
ing that a crossing of different random walks is prohibited. For this purpose 
it is useful to return to Eq. (|8|). It has been shownll3 that the statistics of 
the directed lines with hard-core interaction can be described conveniently 
by replacing the complex field by a field constructed from an algebra of 
nilpotent numbers (i.e., {t^x'")^ = if Z > 1): 

rx^€'" (cT = ±i). 

Each space-time point x is characterized by quantum numbers (a, cr), where 
cr = ±1 corresponds to the real and the imaginary part of This choice 
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implies 2N degrees of freedom at each point x. Therefore, we need also 2N 
variables r/^'"^ at x. Since {vj'^'^} are nilpotent, we have an exclusion principle 
describing the hard-core interaction. Products of ry"'*^ must be commutative 
in order to represent the Bose statistics. Now one can identify an empty site 
with 

and a tunneling process, going from x to x' and connecting the states a and 
a', with 

C a,l a' 2 

j;^Wx,x'Vx Tlx' ' 

where Wx,x' (= Wx~x') describes again the hopping probability as obtained 
from the random-walk expansion. 

The integration over the free Bose field is replaced by a linear mapping 
of the algebra, generated by {^"'"^j, to the complex numbers with 

TT TT TT VT = {^ a A' = A, ix = 1,2,..., N, jx,a = 1,2 _ 

HC Airi,}'^; ^ lo otherwise 

Thus the integral vanishes if the product is incomplete with respect to the 
lattice, the degenerate states a or a. Finally, we must impose periodic 
boundary conditions in the f-direction. It is evident that we can introduce 
analytic functions of the nilpotent field. Therefore, we may write for the 
partition function of the hard-core Bose gas 

f ^ C '2 

x,x' a,a'=l 

The nilpotent field r/"''^ is closely related to a Grassmann field: the algebra, 
generated by {?7"''^}, can be constructed from anticommuting Grassmann 
variablesEl {-i/'"'"", ■0"''^} as 

C''^ = i-iri^ri^x^"- 

Using the usual Grassmann integralEl, the partition function now reads 

Zhc = 

J exp[j2 E i^x,x'6a,a' + J;^Wx,x'm'^P^^'^pir'r^^^^^^ 

(10) 



x,x' a,a'=l 



Here we notice an identity for the diagonal term: 

exp(V'^'iV'x'^C'V°'^) = 1 + c^C'^C'V; 
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The second term on the r.h.s. does not contribute to the partition function, 
since there can only be a factor u^^j.'V'"r^V'"r^ or '^r' ,r'4hf 'ihf fro™ the bosons 
in (|To|). Thus, the diagonal term in Zuc can be replaced by a term that is 
only bilinear in the Grassmann field. The i(7-dependent term in Zhq^ which 
describes the bosons, can also be expressed bilinearly in the Grassmann field 
when we introduce two complex Gaussian fields ("Hubbard- Stratonovich 
transformation"). Two fields are required here, since w is not a positive 
matrix. However, we obtain a positive matrix Vg when we add a positive 
diagonal matrix to w as 

Wx,x' = Wx,x' + {s - s)6x,x' =■ {Vs)x,x' - s6x,x', 

with a sufficiently large s. (Physical results should not depend on s. This 
will be confirmed in the results of the limit — > oo below.) With 

a a 

we obtain the identity 



J exp -^(v9,i;7V)-^(x,x) + (/','^) + (v,P)+^(P,x)+^(x,p)] n^'^^^^^'' 

(11) 

where the scalar product means {f,Vg ^ip) = J2x x' fxVg-x x''^^' ■ "^^^ r.h.s. 
of (11) can be substituted into Zjjc^ ^-^id the bilinear Grassmann term can 



be integrated. This leads to 

/ vrC^ \-P^f f ( rl _i 1 
Zhc = [^) J exp[-N[-{ip,v, ip) + —{x,x) 

- ^'^og[l + {ipx + iXx){'fx+iXx)]\]Yldipxdxx- (12) 

X 

It should be noticed that the field x can also be integrated out in principle. 
However, this is difficult if ^ 1. 

We notice that the expression in the exponent is invariant under a global 
phase transformation of the complex field: (p, x ^ ^^'^'■P^ e*'^x &rid p, X ^ 
e~*'^(^, e~*'^x- Furthermore, the logarithmic term of the action is also sym- 
metric under a local f7(l)-transformation. This means that the depencence 
of the logarithmic term on the local phase of px can be gauged away, and 
this term depends only on the modulus of (fx- We will see in Sect. 5 that 
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the symmetry mider this phase transformation is spontaneously broken. As 
a consequence, there is a one-component Goldstone mode. 

The total density of bosons can be expressed within the new effective 
complex fields 99 and x- From Eq. (|9|) we obtain 



n 



C d 
NMf3dC 



log Zhc 



1 



1 



((^,t;7V)) + -((x,x)) 



+ 



(13) 



where 



(...) 



-m 1 



Zhc 



.. exp 



{-Af[J(V9,i;-V) + ^(x,x) 



^ log[l + {ipx + iXx){Vx + } n d^PxdXx- 



(14) 



3. N ^00: CLASSICAL-FIELD EQUATIONS 



The action in (12) depends on the number of states only through the 
prefactor A^. This suggests that a saddle-point integration can be performed 
for Zhc- For this purpose it is convenient to rescale the fields first: 

We will use the same name for the rescaled fields. Then the saddle-point 
equations read 

dS ^ _i _ ^x + iXx „ 

7^ — = l^i'^s )x,x"fx' - yzm — \(- I ■- ^ = ^ 

d(px ^ C ^ + [fx + lXx)[(Px + IXx) 

^ = Z^K )x',xfx' - 731—7 —■ w- ^ .- ^ = (16) 

9V^x ^ C + [fx + iXx){fx + IXx) 

dS 1_ ijfx + ixx) n n7^ 

C'Xx S C + (Vx + «Xx)(V'a; + ?Xx) 

55 _ 1 ijfx + iXx) _ . . 

OXx s C +[fx+ iXx)[fx + IXx) 

Eqs. (^) - ( |l8|) correspond to a discrete version of the Gross-Pitaevskii (or 
non-linear Schrodinger) equation, often used for the description of a weakly- 
interacting Bose gasQo, a case in which we can assume that fields are small. 
Expansion in terms of the fields yields an equation for ip: 

r 1 1 

In general, it is not possible to solve these equations. Therefore, the assump- 
tion of a uniform solution is useful. This is known as the Thomas-Fermi 
approximationH. 
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3.1. The Thomas-Fermi Approximation 

For a uniform solution we obtain with J2x'i'^s'^)x x' = + and Eqs. 
(0)-® 

— — s s 

There is always a trivial solution <^o = f^o = Xo = Xo = and a non-trivial 
solution with 

= + 

The fugacity has been replaced here by the inverse fugacity z = <^~^. The 
saddle-point contribution from the trivial solution gives for the partition 
function Zq = 1. Any non-trivial solution breaks the symmetry under phase 
transformation for z 7^ 1. The saddle-point equations, however, are invariant 
under the symmetry transformation because the phase factor is not deter- 
mined by them. This leads to massless fluctuations that generate long-range 
correlations. The N 00 limit of the partition function is not affected by 
these fluctuations 

Since the free energy — log Z^c must be minimal, the non-trivial saddle 
point with Zhc = Zi is valid for z < 1 whereas the trivial saddle-point 
solution with Zhc = -^0 is valid for z > 1. This implies for the density of 
bosons of Eq. ( p!^ ) in the limit N ^ 00 



n{z) 



1 — z for z < 1 
for z > 1 



Thus z = 1 is the critical point. According to the definition of z as the 
inverse fugacity of the bosons, z = 1 separates the condensed phase (z < 1) 
from the normal-fiuid phase {z > 1). Corrections of o{N~^) to n(z) will be 
evaluated in Sect. 5.1. 



4. 1/iV-CORRECTIONS: QUASIPARTICLES 

In order to investigate the stability of the saddle-point solutions against 
fluctuations, we shall consider now 

if ^ ipi + 6ip' + i6(p", X^Xi + ^x' + i^x"- 

With Scj) = {5ip' ,6ip" ,6x' ,Sx") we obtain for the partition function in Gaus- 
sian approximation 

Zhc -Z, I exp[-iV((50, g-},5(t>)] J] d6ct>^ (19) 

X 
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with the propagator (Green's function) Qs,v It remains to evaluate the 
latter. For this purpose it is necessary to choose a specific u = e^^ (or 
H = — logti), for instance, 



Ur 



J/2d for r,r' nearest neighbors 

(1-J) for r' = r ■ (20) 

otherwise 



Then the eigenvalues of the matrix w are 

J d 

w = e"^[l-J + -J2cos{kj)] (21) 

with Matsubara frequency lo and d-dimensional wavevector components kj 
with —IT < uj, kj < TT. The eigenvalues are bounded by {wl < 1. It is sufficient 
to consider s = 1 to get a positive matrix vi. Then a simple calculation gives 
the Green's function 



with 

._i 1 fl + k^lAT -uj/2 
1 2V -^^72 l + /tV4r 

for small wavevectors k and frequencies uj. = 0, 1 refers to the two saddle- 
point solutions within the Thomas-Fermi approximation. The new param- 
eter r is proportional to the inverse tunneling rate on the d-dimensional 
lattice: t = d/J. The 2x2 block matrices depend on the fugacity 1/z 

1/1 0\ , /2z-l 
^° = i(o l)' ^^ = ( 1 



5. RESULTS 

The spectrum of the propagator of the field (p is given as (s. App. A) 



<k) = f—^k' + (^?- (22) 

For a weakly-interacting Boge gas of bosons with mass m there is the well- 
known Bogoliubov spectrumEj 



e{k) = ^ nok'^/m + (A;2/2m)2. 
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Since = logC we can identify our parameter r with the mass of weakly- 
interacting bosons. Then Eq. ( p2[ ) agrees with the Bogohubov spectrum 
only for C ~ 1) indicating that the dilute limit of the hard-core Bose gas is 
a weakly-interacting system. The Goldstone mode of Qi is proportional to 

k — uj . 

T 

Thus, the sound velocity is c = \/{C ~ I)/''"- On the other hand, Popovlli 
finds c = a/ fio/m. Therefore, the sound velocity of the N state hard-core 
Bose gas agrees with that of the weakly- interacting bosons only for ~ 1. 



5.1. Total Density of Bosons 



From Eq. (13) we obtain the expression 

n = {I- z)e{l- z) +ni/N + o{N-'^) (23) 

with 

Here Q{x) is the Heaviside step function. This result indicates that there 
are no bosons at z > 1 in the limit N ^ oo, as in the non-interacting Bose 
gas. However, the density increases linearly for z < 1, a regime that is not 
accessible for the non-interacting Bose gas. For a finite number of states 
N there is a non-vanishing density of particles for any value of the inverse 
fugacity z. 



5.2. Density of the Condensate 

The condensate can be studied by considering the correlation function 
of the hard-core Bose field 

1 ^ 

a,a'=l 

On large scales this does not decay in the condensed phase but exhibits 
off-diagonal long range order: 

Cx,x'~no (|r-r'|~oo). 
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0.9 




Fig. 1. Fraction of the condensate no/n as a function of the number of 
degenerate states N per site at fugacity C = 2 (from Eq. (p^)). 



where no is the density of the condensate. Using the notation of Sect. 2 we 
can also express the correlation function of the Bose field as 

Cx,x' = ^^{PxPx')- 

Eventually, in terms of the fields if and x ttie latter reads 

1 ^x + iXx 1 r _, 1 '^x' + iXx' 



I 



'fx' 



N z + {ipx + iXx){fx +%)-' l'''^' N z + {(fx' + iXx')i<fx' + iXx') 

with ip'^ = J2x'{'^s^)x,x'fx'- Using the result of the saddle-point integration 
this expression leads for |r — r'| ^ oo to 

no = (l-^)2(l-z)G(l-z) (24) 

for the density of the condensate. The prefactor (1 — 1/N)'^ indicates a 
suppression of the condensate due to increasing interaction with a decreasing 
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number of particles states N at each lattice site. In Eq. (|24| ) terms of higher 
order in are neglected. It is expected that they lead to additional 
terms such that no does not vanish at = 1. Combining ( p3|) and the 
condensate fraction reads 

which is plotted in Fig.l. 



6. CONCLUSIONS 

A hard-core Bose gas on a lattice has been used to study the effect 
of interaction in an optical lattice. Introducing A'^ degenerate states at each 
lattice site and a hard-core interaction only between the same state, the limit 
N ^ oo can be solved by a saddle-point integration. This limit exhibits a 
transition from a normal state to a Bose-Einstein condensate if the fugacity 
C = 1. The limit A/^ ^ oo is very special, since the total density of bosons 
vanishes for C < 1- However, 1/A^-corrections lead to a non- vanishing density 
of bosons for any value of (. On the other hand, 1/A^-corrections indicate a 
depletion of the condensate {( > 1). These results agree with the observation 
that the condensate fraction decreases with an increasing coupling constant 
in the Bose-Hubbard mo del. Further studies are neccessary in order to find 
the effect of very strong interaction when ~ 1. 
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APPENDIX A: EFFECTIVE PROPAGATOR FOR THE FIELD 

f 

Starting from the partition function Zhc of Eq. the field x can 

be integrated out. The result is an effective propagator of the field (p 

Qv 

where we have 



0\ \({2z-\flz 

ij' ^^ = 2 1 1 
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Then a straightforward calculation gives the following expressions 



and 



1 I ({z-l)/{z + l) + e/AT -u/2 

-u/2 (z - l)/(z + 1) + A;V4t 



1 I ( {I- z)/z + k'^/AT -u/2 



-u/2 k^l^T^ 

Prom Q^^ we obtain for the spectrum of the condensed phase (i.e. z < 1) 
the expression 
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